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OBJECTIVES 

Oui  overall  goal  )s  to  develop,  implement  and  transfer 
new  nometical  methods  for  solving  the  moving 
ixjtdfacc  problems  of  solidihcation  smd  crystal  growth*  We 
have  sevOTl  speci&c  objectives; 

•  Combine  fest  algorithms,  level  set  techniques,  adap¬ 
tive  refinement  and  data-  stmctnres  to  develop  and 
implement  accurate,  efficient  and  general  new  meth¬ 
ods  for  moving  sharp  interfeces. 

•  Develop  a  modular  level  set  code  for  transfer  to  other 
researchers,  labs,  and  industry. 

•  Combine  spectral  methods  and  stiff  ODE  solvers  to 
produce  fast  spectrally-accurate  solvers  for  phase  field 
models  of  solidification  and  related  problems. 

•  Validate  phase  field  models  by  numerical  experiment 

and  identify  parameter  choices  appropriate  to  indus¬ 
trial  settii^. 

STATUS  OF  EFFORT 

We  have  attained  our  first  and  second  objectives  with 
the  modular  semi-Lagrangian  methods  for  moving  inter¬ 
faces  reported  in  Publications  [P2-PS].  Our  third  objec¬ 


tive  was  completed  for  the  case  of  periodic  boundary  condi¬ 
tions  in  (Strain,  1995b),  whici  is  sufficient  for  basic  validar 
tion  and  parameter  identification*  Further  work  on  general 
boundary  conditions  and  general  domains  is  in  progress. 
Additionally,  we  present  effident  and  accurate  new  vortex 
methods  for  convection  in  the  mdt  in  Publication  [Pi]. 

ACCOMPLISHMENTS  AND  NEW  FINDINGS 

We  now  summarize  our  progress  in  vortex  methods  and 
moving  interfaces  in  more  detril. 

Vortex  methods 

We  have  studied  vortex  methods  for  incompressible  2D 
flows  both  as  the  simplest  examples  of  level-set-based  mov¬ 
ing  boundary  techniques  and  as  solvers  for  the  convection 
in  the  mdt  around  a  growing  crystal.  Our  most  recent 
work  on  vortex  methods  is  concisely  s\irveyed  in  (Strain, 
1996b). 

♦  ia  1992-4,  we  developed  a  new  second-order  accu¬ 
rate  ‘triangulated’^  vortex  method  which  solved  the  long¬ 
standing  problem  of  long-time  inaccuracy  in  traditional 
vortex  methods  (Russo  and  Strain,  1994). 

♦  In  1994-5,  we  carried  this  line  of  research  further. 
First,  we  developed  an  important  tool  for  constructing 
quadrature  formulas  for  singular  integral  operators  in  sev¬ 
eral  dimensions  (Strain,  199Sa). 
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•  Then  in  1996,  we  produced  new  quadrature-based  2D 
vortex  methods  of  arbitrary  order  with  excellent  long-time 
accuracy  properties  (Strain,  1996a), 

•  In  1996,  we  also  developed  new  methods  based  on 
the  smooth  quadrature  roles  developed  in  (Strain,  1995a), 
Publication  pi]  presents  an  accurate  and  efficient  new  ap- 
proadi  to  vortex  methods  for  modeling  incompressible  in- 
visdd  fluid  flow  in  two  space  dimensions. 

Abstraa 

Since  thdr  introduction  in  the  1970’s,  vortex  methods 
have  been  widely  generalized  and  ^plied  to  compute  many 
complex  physical  flows.  They  model  sudi  flows  by  moving 
discrete  particles  with  a  smoothed  ^proximate  fluid  veloc¬ 
ity.  Thus  thdr  convergence  analysis  must  balance  smooth¬ 
ing  errors  against  velodly  approximation.  Three  long¬ 
standing  theoretical  problems  have  hampered  the  analysis 
of  vortex  methods: 

o  Long-time  computations  lose  high-order  accuracy  via 
the  Perlman  effect:  velodty  ^proximatioa  errors 
blow  up  when  the  flow  map  twists. 

o  Proofe  of  order-m  convergence  require  2m  deriva¬ 
tives,  imidi  more  than  standard  numerical  methods. 

o  Parameter  balances  a{>propriate  for  convergence  pre¬ 
clude  effident  velocity  evaluation  by  the  fast  multi- 
p(^  method. 

Reformulation,  quadrature  and  error  analysis  are  used 
in  this  paper  to  solve  these  problems.  The  Perlman  effect 
is  eliminated  hy  combining  a  freo-Lagrangian  formulation 
with  fast  adaptive  quadratiire  rules.  Thus  long-time  high- 
order  accuracy  is  obtained  by  constructing  a  new  set  of 
particle  weights  at  each  step.  A  new  error  analysis  of  ve¬ 
locity  approximation  suggests  a  balance  of  parameters  and 
choice  of  smoothing  functions  whi^  requires  only  m  -f*  1 
derivatives  for  essentially  order-m  convergence,  halving  the 
differentiability  required  of  the  flow.  This  balance  can  be 
combined  with  &st  sommation  schemes  to  yield  almost 
optimal  efficiency.  Numerical  results  verity  the  error  anal¬ 
ysis  and  erfdbit  excellent  accuracy  and  efficiency,  giving 
this  approach—whidi  generalizes  widely — practical  value 
as  well  as  theoretical  importance. 

Summary 

Vortex  methods  involve  several  components:  velodty 
evaluation,  vortex  motion,  diffusion,  boundary  conditions 
and  regridding.  In  Publication  [Pi],  we  improve  the  speed, 
accuracy  and  robustness  of  the  velocity  evaluation.  We 
eliminate  the  flow  map,  improve  the  quadrature  used  for 


the  Biot-Savait  law,  and  analyze  the  error  in  a  nonstan¬ 
dard  way,  requiring  less  differentiability  of  the  flow  and 
obtaining  effident  new  parameter  balances.  We  employ 
standard  tedmiques  for  the  vortex  motion  and  consider  in- 
visdd  free-space  flow  to  eliminate  diffusion  and  boundary 
conditions.  Our  approach  combines  naturally  with  regnd- 
diTig  and  fest  summation  methods. 

First,  we  redew  Lagrangian  vortex  methods.  These 
move  the  nodes  of  a  fixed  quadrature  rule  with  the  com¬ 
puted  fluid  velodty,  preserving  the  weights  of  the  rule 
by  incompressibility-  This  procedure  loses  accuracy  when 
the  flow  becomes  disorganized,  motivating  many  regnd- 
ding  techniques.  Even  before  the  flow  becomes  disorga- 
nized,  however,  obtaining  high-order  accuracy  with  a 
quadrature  rule  requires  smoothing  of  the  singular  Biot- 
Savart  kernel.  Smoothing  high-^rder  accuracy  for 
short  times  but  slows  down  fast  velodty  evaluation  teAr 
niques  and  halves  the  order  of  accuracy  relative  to  the  dif¬ 
ferentiability  of  the  flow.  We  review  two  free-Lagrangjan 
vortex  methods,  the  triangulated  vortex  method  we  devel- 
Of>ed  in  (Russo  and  Strain,  1994)  and  the  quadrature-based 
method  of  (Strain,  1996a).  Triangulated  vortex  methods 
are  robust,  practical  and  effident  but  limited  to  second- 
order  accuracy.  Quadrature-based  methods  compute  adajy 
tive  quadratures  tailored  to  the  Biot-Savait  kernel  at  each 
time  step,  yidding  long-time  hi^h-order  accuracy  at  atymp” 
totically  optimal  cost. 

Publication  [PI]  devdpps  a  firee-Lagran^an  method  whi  A 
couples  kernel  smoothing  wiA  adaptive  quadrature  rules 
not  tailored  to  Ae  Biot-Savait  fcem^  producing  long-time 
hi^-order  accuracy.  The  asymptotic  slowdown  produced 
by  kemd  smoothing  is  almost  duninated  by  a  careful  chwee 
of  smoothing  functions  and  parameters,  based  on  a  new 
error  analysis  of  the  velocity  evaluation.  This  analysis  re¬ 
quires  about  half  as  many  derivatives  of  Ae  solution  as  the 
standard  approaA. 

The  structure  of  our  mcAod  is  standard:  At  eaA  time 
step,  the  smoothed  velodty  is  evalxiated  once  and  Ae  vor¬ 
tices  are  moved  wiA  an  explidt  multistep  meAod.  The 
velodty  evaluation  is  nonstandard;  First,  a  data  structure 
groups  Ae  N  vortices  into  cells  convenient  for  integration. 
Then  a  global  order-q  quadrature  rule  is  built.  Finally,  Ae 
fest  multipole  method  is  used  wi A  this  rule  to  evaluate 
the  smoothed  velodty  field. 

The  error  Is  measured  for  standard  test  problems  and 
otir  theoretical  predictions  are  fiilly  verified.  Then  more 
comjJex  flows  are  computed- 

Moving  inteifaces 

Our  work  on  moving  interface  problems  in  matenals 
science  combines  fast  PDE  solvers  suA  as  boundary  int^ 
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gral  methods  with  fast  geometric  algorithms  a&d  level  set 
techniques  to  build  efiFective  new  numerical  methods, 

•  In  19S9,  we  developed  a  powerflil  implicit  boundary 
integral  method  for  computing  periodic  dendrite  forma¬ 
tion  in  a  Stefan-t>pe  model  of  unstable  solidification.  We 
used  a  spectral  discretization  of  a  new  arclength-preservlDg 
curve  motion  algorithm  and  an  unoonditionally-stable  im¬ 
plicit  time  stepping  scheme,  solved  by  a  modified  New¬ 
ton  method.  method  was  expensive  but  accurate  and 
robust;  higherorder  accuxa^r  was  obtained  by  Richard¬ 
son  extrapolation.  Numerical  results  agreed  closely  with 

stability  theory  and  predicted  correct  tip-splitting 
phenomena  in  the  nonlinear  regime.  The  method  and  nu¬ 
merical  results  are  presented  in  (Strain,  19S9). 

•  In  1990,  we  devdoped  fast  algorithms  for  evaluat¬ 
ing  heat  potentials  (Greengard  and  Strain,  1990)  which 
speeded  up  our  boundary  integral  method  (and  many  other 
^culadons)  by  several  orders  of  magnitude.  This  dis¬ 
covery  led  to  the  fost  Gauss  transform  which  has  proven 
widely  explicable  (Greengard  and  Strain,  1991;  Strain, 

1991) , 

•  In  1990-2,  we  developed  a  level  set/boundary  integral 
method  for  dendritic  solidification  (Sethian  and  Strain, 

1992) .  This  combined  the  boundary  integral  method  of 
(Strain,  1989)  with  fast  algorithms  from  (Greengard  and 
Strain,  1990)  and  the  level  set  method  of  (OAer  and  Sethian, 
1988).  The  levd  set  method  bandied  topolo^cal  changes 
e&ctively  while  fast  boundary  integral  techniques  ensured 
accuracy  and  efficiency  in  the  vdocity  evaluation, 

•  In  1992-6,  we  developed  efficient  and  accurate  new 
vortex  methods  for  modeling  convection  in  the  melt  (Russo 
and  Strain,  1994;  Strain,  1996a),  together  with  new  error 
analyses  [Pi]  and  quadrature  rul©  (Strain,  1995a)  for  gen¬ 
eral  integral  equations.  Much  of  this  work  is  summarized 
above  and  in  (Strain,  1996b). 

•  Our  major  achievements  in  1996-99  have  been  the 
devdopment  and  implementation  of  hi^y  effective  new 
numerical  methods  for  general  moving  interface  problems. 
Abstracts  of  Pubfications  [P2-P5],  then  an  overview  of 
these  methods,  follow. 

Semi-Lagrangjao  methods  for  level  set  equations  (P2] 

A  new  numerical  method  for  solving  geometric  moving 
interface  problems  is  presented.  The  method  combines  a 
levd  set  approach  and  a  semi-Lagrangian  time  stepping 
scheme  which  is  explicit  yet  unconditionally  stable.  The 
combination  decouples  each  mesh  point  from  the  others 
and  the  time  step  from  the  CFL  condition,  permitting  the 
construction  of  methods  which  are  efficient,  ad^tive  and 
modular. 

Analysis  of  a  linear  one-dimendonal  modd  problem 


suggests  a  surprising  convergence  criterion  which  is  sup¬ 
ported  by  heuristic  arguments  and  confirmed  by  an  exten¬ 
sive  collection  of  two-dimensional  numerical  results.  The 
new  method  computes  correct  viscosity  solutions  to  prob¬ 
lems  involving  geometry,  anisotropy,  curvature  and  com¬ 
plex  topological  events. 

Fast  tree-based  redistandng  for  levd  set  computations  [P3] 

Level  set  methods  for  moving  interface  problems  re¬ 
quire  efficient  techniques  for  transforming  an  interface  to 
a  globally  defined  function  whose  zero  set  is  the  interface, 
sudi  as  the  signed  distance  to  the  interface. 

This  paper  presents  an  efficient  algorithm  for  this  ‘‘re- 
distandn^  problem.  The  algorithm  uses  trees  and  fast 
Ddaunay  triangulation  to  compute  a  global  approximate 
Signed  distance  function.  An  adaptive  tree  mesh  is  built 
to  resolve  the  interface,  and  the  vertex  distances  are  eval¬ 
uated  exactly  with  an  efficient  seardi  strategy,  providing 
both  continuous  and  discontinuous  interpolants.  Given  a 
polygonal  interface  with  N  elements,  the  algorithm  runs 
in  O^Nlo^N)  sp&ce  and  time,  and  numerical  results  show 
it  is  highly  efficient  in  practice. 

Tree  rr^ethods  for  moving  interfaces  [P4] 

A  fast  adaptive  numerical  method  for  solving  moving 
interface  problems  is  presented.  The  method  combines  a 
level  set  approadi  with  frequent  redistandng  and  a  sezni- 
Lagrangjan  time  stepping  scheme  which  is  explidt  yet  un¬ 
conditionally  stable.  An  adaptive  tree  me^  is  used  to 
concentrate  computational  effort  on  the  interface,  so  the 
method  moves  an  interface  with  N  degrees  of  freedom  in 
©(NlogJ'f)  work  per  time  step,  Effidency  is  increased 
by  taking  large  time  steps  even  for  parabolic  curvature 
flows.  The  method  computes  accurate  viscosity  solutions 
to  a  wide  variety  of  difficult  moving  mtei&oe  problems  in¬ 
volving  merghg,  anisotropy,  faceting  and  curvature. 

A  fast  modular  semi-Lagrangian  method  for  moving  interfaces 

A  fast  modular  numerical  method  for  solving  general 
moving  interfooe  problems  is  presented.  It  simplifies  code 
devdopzneat  by  providing  a  blade-box  solver  whidi  moves 
a  ^ven  interface  one  $t^  with  given  normal  velod^.  The 
method  combines  an  effidently  redistanced  level  set  ap¬ 
proach,  a  problem-independent  velodty  extension,  and  a 
second-order  semi-Lagrangian  time  stepping  scheme  wluch 
reduces  numerical  error  by  exact  evaluation  of  the  rigned 
distance  function. 

Adaptive  quadtree  meshes  are  used  to  concentrate  com¬ 
putational  effort  on  the  interface,  so  the  method  moves  an 
N'-dement  interface  in  0{N  log  N)  work  per  time  step. 
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Efficiency  is  increased  by  taking  large  time  steps  even  for 
parabolic  curvature  flows.  Numerical  res\ilts  show  that  the 
method  compute  accxirate  viscosity  solutions  to  a  wide 
variety  of  difficult  geometric  moving  interface  problems  in¬ 
volving  merging,  anisotropy,  faceting,  nonlocality  and  cur¬ 
vature. 

Definitions  and  examples 

A  moving  interface  is  a  collection  r(t)  of  nonintersect- 
ing  oriented  closed  curves  (in  the  plane)  or  surfaces  (in 
space).  A  suffidently  smooth  moving  interface  has 

©  An  outward  unit  normal  vector 

©  A  signed  curvature  C{x,t),  and 

o  A  normal  velocity  function  V(x,t). 

A  moving  interface  problem  consists  of  a  specification  of 
V (z,  t)  as  a  function  of  r(f),  its  history  and  geometry,  and 
any  other  fields  or  variables  which,  may  be  present.  For 
example,  we  may  have  passive  transport 

under  a  ^ven  velocity  field  F,  or  curvature-  and  angle- 
dependent  normal  velocity 

V(x,  t)  =  H  +  £COs{*(5  +  5o))  + 


Level  set  methods 

The  level  set  method  (Osher  and  Sethian,  1988)  turns 
around  the  last  formula  for  V  and  regards  it  as  a  PDE  for 

V^-V|lVs^lI=0. 

This  requires  us  to  extend  V  off  r(f )  to  a  function  defined 
everywhere.  This  “V'  extension’*  is  one  of  the  ingredi¬ 
ents  of  the  level  set  tedmique.  The  V  extendon  should  be 
smooth  and  should  agree  with  V  on  r(f),  but  is  otherwise 
arbitrary.  The  numerical  properties  of  the  method  can  be 
stron^y  affected  by  the  choice  of  V  extension. 

The  standard  levd  set  method  involves  four  steps. 

©  Choose  a  V  extension. 

©  Evaluate  ^  on  a  uniform  mesh. 

©  Advance  (p  with  a  uniform  mesh  solver  borrowed 
from  hyperbolic  conservation  laws. 

o  Find  r(t)  by  contouring  tp  when  dedred. 

A  fifth  step,  redistandng,  can  be  performed  occasionally 
to  enhance  robustness; 


where  cos  5  =  N* •£,  or  we  may  have  Stefan-type  conditions 
such  as 

n...)  =  [§ 

where  u(x,  t)  solves  the  heat  equarion  Ut  =  Au  off  r(t)  and 
tt  =  eC  on  F(t),  as  in  £ar  dendritic  solidification  of 

a  pure  material  such  as  sucdnonitrile. 

Level  set  equations 

Any  moving  interface  problem  can  be  reformulated  as 
a  PDE  for  a  fimction  <p  whose  zero  act  is  r(f): 

r(f)  =  {x\ip{x,t)  =  0}, 

for  example  the  dgned  distance  function 

with  the  +  sign  diosen  in  the  interior  of  r(t).  The  normal, 
curvature  and  velocity  are  given  by 

N=^ 

livin’ 


o  R^lace  <p  by  the  dgned  distance  to  its  zero  set  r(t). 

On  a  uniform  mesh,  redistandng  can  be  extremdy  costly — 
as  costly  as  the  rest  of  the  calculation  put  together. 

Advantages  and  disadvantages 

The  main  advantage  of  the  levd  set  method  is  topology 
icol  rohistn^ss:  meri^ng  and  splitting  pieces  of  interface  is 
handled  automatically  by  the  levd  set  equation.  Potential 
disadvantages  include  the  expense  of  going  up  a  diznen- 
don,  the  grid-dependence  of  the  solution,  and  the  diffi¬ 
culty  of  flrtdipg  a  problem-dependent  V  extendon.  Since 
the  method  is  es^lidt ,  tiny  time  steps  are  required  for  some 
problems.  Stefan-type  problcns,  for  example,  require  the 
time  step  k  to  satisfy 

k  <  0{h^) 


on  a  mesh  of  size  A. 

Our  new  method  overcomes  these  disadvantages  by  com¬ 
bining  the  levd  set  approach  with  an  adaptive  tree  mesh 
for  efficiency,  a  general  problem-independent  V  extension, 
and  an  explicit  unconditionally  stable  time  stepping  method. 
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Time  stepping 

Our  time  stepping  method  is  based  on  the  backward 
characteristic  method  of  Courant,  Isaacson  and  Rees  (1953). 
The  “CIR^  method  sol^res  the  hyp^bolic  equation 

9  -  F(x,  i)  •  V<p  =  0 

by  the  following  algorithm:  At  each  z  in  the  grid, 

o  Evaluate  the  velocity  F(z,tn)- 

o  Move  X  back  with  velocity  F(^,  tn)  to  s  =  X’hkF{x,  t^). 

o  Interpolate  <p{x,tn)  to  the  point  s. 

o  Set  equal  to  the  inteipolated  value. 

The  CIR  method  is  first-order  accurate  and  not  in  con- 
servation  form,  but  has  the  unusual  property  of  being  ex¬ 
plicit  yet  uncondidonallj  stable.  It  converges  fi>r  model 
problems  as  long  as 

k>0{h),  . 

an  accuracy  condition  which  is  the  reverse  of  the  usual 
CFL  stability  condition. 

We  have  developed  a  second-border  “semi-Lagrangian" 
method  which  combines  a  C2R  predictor 

^(®)  =  0  =  v»(« + *■?(*>  t).  0  (I) 

vith  a  tiapezoidal  corrector 

^(*)  =  <fiix  +  t)  + t  +  (2) 

Here  F  is  built  &om  the  interfacia]  velocity  V  of  the  zero 
set  GaTTzma  of  the  predicted  solution  ^  at  time  f -I-  This 
predictor^-corrector  pair  is  second-order  accurate  in  time, 
explicit,  and  unconditionally  stable;  eadi  time  step  re¬ 
quires  two  velocity  evaluations.  Since  our  advection  veloc¬ 
ity  F(x,  t)  extends  the  user-specified  velocity  functional  V 
on  the  zero  set  r(t)  of  t),  each  semi-Lagran^an 

time  step  requires  several  complex  ^obal  operations.  Start¬ 
ing  with  an  interface  F(t)>  our  method  carries  out  the  fol¬ 
lowing  steps  to  produce  the  new  interface  r(f  +  fc): 

o  Evaluate  the  signed  distance  <p  from  the  interface 

r(t). 

6  Evaluate  the  interfedal  velocity  V  of  T{t)  by  a  user- 
supplied  module. 

o  Extend  V  to  a  ^obal  advection  velocity  jF. 

o  Advance  via  F  to  the  predicted  CIR  solution  tp 
defined  by  Eq.  (1). 


o  Contour  ^  to  get  the  predicted  interface  f . 

o  Evaluate  the  predicted  interfacial  velocity  V  of  f. 

o  Extend  V  to  a  global  advection  veloaty  P. 

o  Advance  via  F  and  F  to  the  corrected  solution  tp 
defined  by  Eq.  (2). 

o  Contour  ^  to  get  r(f  -r  i). 

Each  of  these  steps  can  be  efficiently  implemented  with  an 
adaptive  quadtree  mesh. 

Semi-Lagrangian  methods  for  level  set  equations 

We  first  tested  the  CIR  scheme  in  Publication  [P2]  by 
solving  level  set  equations  with  a  fixed  uniform  messh  and 
ENO  differendng,  to  obtain  simplicity  and  hi^-order  a& 
curacy.  We  extended  the  velocity  for  passive  transport 
and  geometric  moving  interfeces  by  evaluating  the  natu¬ 
ral  formulas,  with  frequent  redistancing,  smoothing  and 
truncation  sufficient  to  satisfy  the  CFL  conffition  even  for 
parabolic  problems  like  motion  by  curvature.  A  surprising 
result  of  tlds  work  is  that  the  CFL  condition  can  be  satis¬ 
fied  by  nonlocal  velocity  evaluation  rather  than  restricting 
the  time  step.  Numerical  results  verify  the  accuracy  and 
stability  of  the  method. 

Fast  tree-based  redistancing  for  level  set  computatiorts 
As  an  intermediate  step  between  the  uniform  mesh 
method  and  our  ultimate  goal.  Publication  [P3]  develops 
an  effident  new  redistanong  technique  with  the  aid  a 
new  data  structure  called  the  distance  tree.  This  struc¬ 
ture  is  a  quadtree  whose  cells  know  their  distance  to  r(t) 
and  nearby  elements  of  r(t),  allowing  for  effident  approxi¬ 
mate  redistandngof  Figure  1  shows  a  pentagonal  curve 
with  the  corresponding  &Ievd  distance  tree  and  signed 
distance  function.  The  distance  tree  is  built  effidently 
by  a  three-step  recursiye  seardi  procedtire  with  guaran¬ 
ty  correctness.  Ddauziay  triaagulation  is  used  to  fix  the 
sign  of  the  signed  distance  function  The  technique  is 
asymptotically  optimal,  requiring  only  (OiV' log  N)  work 
to  redistance  an  interface  with  N  elements,  and  extremely 
in  practice:  computations  of  moderate  complexity  are 
speeded  up  400  times,  while  redistanong  the  CIR  calcula¬ 
tion  on  a  uniform  mesh  costs  considerably  less  than  morixig 
tp  one  step. 

Tree  methods  for  moving  interfaces 

Our  work  with  distance  trees  leads  in  to  the  adaptive 
tree-based  level  set  method  of  Publication  [P4].  Here  we 
combine  backward  characteristics  with  an  adaptive  tree 
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mesh  to  build  a  method  which  is  not  only  accurate  and  ro¬ 
bust  but  also  optimally  efficient:  An  interface  T{t)  with  N 
degrees  of  freedom  costs  only  0{N  logN)  to  move  one  step. 
Efficiency  is  further  enhanc^  by  the  semi-Lagrangiaa  time 
stepping  scheme:  Large  time  steps  can  be  taken  even  though 
the  tree  mesh  contains  tiny  cells,  because  the  backward 
diaracteristic  method  is  unconditionally  stable. 

The  tree  mesh  is  refined  with  a  new  functional  ap¬ 
proach:  Given  a  level  set  function  ^  build  a 

tree  at  time  -4-  i  by  recursive  evaluation  of 

tn+i)  =  tn)  at  projected  points  $  =  x-i^kFix,  t«). 
The  criterion  for  splitting  a  tree  cell  is  Ampler  the  values  of 
fp{x^  tn+i)  on  the  cell  are  smaller  than  the  size  of  the  cell. 
This  contrasts  favorably  with  standard  mesh  refinement 
algorithms,  which  tend  to  be  based  on  error  estimates  for 
the  computed  solution. 

In  our  numerical  experiments,  the  velocity  is  evaluated 
either  directly  cm  the  tree  or  by  transferring  9  to  a  uniform 
grid  and  empl<tying  the  velocity  evaluation  tedmique  of 
Publication  [P2].  The  sample  calculations  shown  in  Rgure 
2  clearly  demonstrate  the  robustness  of  the  tree  method, 
even  when  computing  complex  merges  between  interfaces 
with  comers-  Grid  ^ects,  which  often  hamper  level  set 
methods,  are  clearly  absent. 

A  fast  modular  semi-Lagrangian  method  for  moving  inter¬ 
faces 

A  completely  modular  method  for  moving  interfaces 
is  developed  and  validated  in  Publication  [P5].  It  com¬ 
bines  the  ideas  above  with  second-order  semi-lagrangian 
time  stepping,  efficient  exact  qiiadtree-based  redistandng, 
a  modular  problem-independent  velocity  extenson,  and 
exact  ^interpolation  in  the  semi— Lagran^an  method.  Our 
vdocity  extension  tedmique  evaluates  the  nearest-point 
extension  on  a  distance  tree,  builds  a  contimious  inter- 
polant,  and  satisfies  a  Tnaximiim  prindple. 

The  modular  method  resolves  and  moves  complex  in¬ 
terfaces  at  optimal  cost  with  time  steps  unconstrained  by 
numerical  stability.  It  is  a  “black-box”  method  for  moving 
interfaces,  which  accepts  the  interface  and  its  velodty  at 
time  t  and  returns  the  evolved  interface  one  time  step  later. 
Such  methods  simplify  moving  interfaces,  because  the  nu¬ 
merics  are  independent  of  the  phyrical  problem  driving  the 
interfadal  motion. 

Numerical  results  show  that  our  method  converges  to 
correct  viscosity  solutions  even  for  difficult  moving  inter¬ 
face  problems  involving  mer^ng,  faceting,  transport,  non- 
locafity  and  anisotropic  curvature-dependent  geometry.  Large 
time  steps  can  be  taken  even  for  parabolic  problems  such  as 
curvature  flows.  Selected  computational  results  are  shown 
in  Figure  3. 
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Figure  1:  Six-levd  distance  tree  (a)  and  levd  set  signed- 
distance  function  on  tree  (b). 


Figure  2:  Sample  solutions  of  moving  interface  problems: 
(a)  initially  circular  bubbles  after  transport  in  a  shear¬ 
ing  flow,  (b)  merging  of  complex  interfiaces  with  unit  nor¬ 
mal  vdocity*  and  (c)  crystalline  facets  developing  under  a 
threefold  anisotropic  curvature-dependent  normal  velodty. 
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Figure  3:  Selected  computational  results:  (a-b)  passive  tranq>ort,  (c)  unit  normal  velociQr,  (d-«)  shrinking  and  merging 
into  Wolff  shapes,  (f)  curvature  flow,  (g-i)  spiral  and  (j-1)  trefoils  under  volume-preserving  curvature  flow. 
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